Efficient quantification of non-Gaussian spin distributions 
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We study theoretically and experimentally the quantification of non-Gaussian distributions via 
non-destructive measurements. Using the theory of cumulants, their unbiased estimators, and the 
uncertainties of these estimators, we describe a quantification which is simultaneously efficient, unbi- 
ased by measurement noise, and suitable for hypothesis tests, e.g., to detect non-classical states. The 
theory is applied to cold 87 Rb spin ensembles prepared in non-gaussian states by optical pumping 
and measured by non-destructive Faraday rotation probing. We find an optimal use of measurement 
resources under realistic conditions, e.g., in atomic ensemble quantum memories. 
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Introduction - Non-Gaussian states are an essential re- 
quirement for universal quantum computation [l|, Q and 
several quantum communication tasks with continuous 
variables, including improving the fidelity of quantum 
teleportation Q and entanglement distillation 

HI!- Op- 
tical non-Gaussian states have been demonstrated [6|- 
[Toj and proposals in atomic systems (llTflij are being 
actively pursued. In photonic s yste ms, histograms 
and state tomography 0, @, foL llCf have been used 
show non-Gaussianity, but require a large number of mea- 
surements. For material systems with longer time-scales 
these approaches may be prohibitively expensive. Here 
we demonstrate the use of cumulants, global measures of 
distribution shape, to show non-Gaussianity in an atomic 
spin ensemble. Cumulants can be used to show non- 
classicality flil- 18|. can be estimated with few measure- 
ments and have known uncertainties, a critical require- 
ment for proofs of non-classicality. 

Approach - Quantification or testing of distributions has 
features not encountered in quantification of observables. 
For example, experimental measurement noise appears as 
a distortion of the distribution that cannot be "averaged 
away" by additional measurements. As will be discussed 
later, the theory of cumulants naturally handles this sit- 
uation. We focus on the fourth-order cumulant K4, the 
lowest-order indicator of non-Gaussianity in symmetric 
distributions such as Fock and "Schrodinger kitten" 
states @, EH ■ We study theoretically and experimentally 
the noise properties of Fisher's unbiased estimator of K4, 
i.e., the fourth "k-statistic" ki, and find optimal measure- 
ment conditions. Because K4 is related to the negativity 
of the Wigner function this estimation is of direct 
relevance to detection of non-classical states. We employ 
quantum non-demolition measurement, a key technique 
for generation and measurement of non-classical states 
in atomic spin ensembles 2^, 21 and nano- mechanical 
oscillators [22|. 

Moments, cumulants and estimators - A continuous 
random variable X with probability distribution func- 



is completely characterized by its mo- 
= J X k P(X)dX or cumulants K n — \i n — 
Mn-fe K fci where (J?) is the binomial coeffi- 



tion P(X) 
ments /ifc = 

Z^fc=l \, 
cient. 

Since Gaussian distributions have n n>2 = 0, estima- 
tion of K4, (or K3 for non-symmetric distributions), is a 
natural test for non-Gaussianity. In an experiment, a fi- 
nite sample {X\ . . .Xn} from P is used to estimate the 
k's. Fisher's unbiased estimators, known as "k-statistics" 
k n , give the correct expectation values (k n ) = K n for fi- 
nite N [13. Defining S n = ^ X f we have: 

k 3 = (2Sf - 3NS!S 2 + N 2 S 3 )/N {2 



(1) 



k 4 = (-6Sf + 12NSfS 2 - 3N(N - l)S% 



-4iV(iV-l)S , iS , 3 + A 2 (iV + l)5'4)/A ( 3 ) (2) 

where N {m} = N(N - 1) ... (N - m) . 

We need the uncertainty in the cumulant estima- 
tion to test for non-Gaussianity, or to compare non- 
Gaussianity between distributions. For hypothesis test- 
ing and maximum-likelihood approaches, we need the 
variances of k 3 ,k 4 for a given P. These are found by 
combinatorial methods and given in reference [23| : 

var(fc 3 ) = kq/N + 9JV(k 2 k 4 + K 2 )/iV (1) + QN 2 nl/N (2) (3) 
var(fc 4 ) = k$/N + 2N(8k 6 k 2 + 24k 5 k 3 + 17k 2 )/A (1 ) 

+72A 2 (k 4 K2 + 2k 2 k 2 )/A (2) 

+24N 2 (N + 1)4/N (3) . (4) 

It is also possible to estimate the uncertainty in k 4 from 
data {X} using estimators of higher order cumulants [23] ]. 
The efficiency of cumulant estimation is illustrated in Fig. 

HI 

Measurement noise - When the measured signal is Z — 
X + Y, where X is the true value and Y is uncorre- 
cted noise, the measured distribution is the convolu- 
tion P{Z) = P(X) ® P(Y). The effect of this dis- 
tortion on cumulants is the following: for independent 
variables, cumulants accumulate (i.e., add) so that 



2 




FIG. 1. (color online) Simulated estimator ki as a function 
of sample size N. a) (insets) black curves show quadrature 
distributions of states p = (1 — p) 0) (0| (1|, scaled to 

unit variance, and six N = 1000 histograms (offset for clarity) 
for p = (green), 1/3 (brown), 1/2 (blue) and 2/3 (red), b) 
Ten realizations of ki versus N drawn from each of the four 
distributions. Shaded regions show K4 ± ^/var(Ai4), from Eqs 
©, ©. With N = 1000, fc 4 distinguishes p = 1/2 (blue) 
from p = (green, Gaussian) with > 7<j significance, even 
though the histograms look similar "to the eye." 



„(*) 



,X) 



where , k^ indicate k„ 



k n for 



distribution P(Q). The extremely important case of un- 
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and 

(X) 



correlated, zero-mean Gaussian noise, k 2 
other cumulants zero, is thus very simple: 
except for k 2 = K 2 + °y- Critically, added Gaussian 
noise does not alter the observed K3, K4. 
Experimental system and state preparation - We test this 
approach by estimating non-Gaussian spin distributions 
in an atomic ensemble, similar to ensemble systems being 
developed for quantum networking with non-Gaussian 
states [24|]. The collective spin component F z is mea- 
sured by Faraday rotation using optical pulses. The 
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GN L F z /2, 
TVl is the number of 



detected Stokes operator is S y 
where G is a coupling constant 
photons, and S y is the input Stokes operator, which 
contributes quantum noise. In the above formulation 



X = F z , Y = 2S { y in) /(GN L ) and Z = 2S^'/{GN L ) 

The experimental system is described in detail in ref- 
erences 
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26|. An ensemble of ~ 10 6 87 Rb atoms is 
trapped in an elongated dipole trap made from a weakly 
focused 1030 nm beam and cooled to 25 uK. A non- 
destructive measurement of the atomic state is made us- 
ing pulses of linearly polarized light detuned 800 MHz to 
the red of the F = 1 — > F' — transition of the D2 line 
and sent through the atoms in a beam matched to the 
transverse cloud size. The pulses are of 1 us duration, 
contain 3.7 x 10 6 photons on average, and are spaced 
by 10 us to allow individual detection. The 240:1 aspect 
ratio of the atomic cloud creates a strong paramagnetic 
Faraday interaction 

G « 6 x 10~ 8 rad/spin. After interaction with the 
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Baseline Displace and Count 
acquisition measure atoms 

FIG. 2. Experimental sequence: The experimental sequence 
divides into distinct tasks. Baseline acquisition: prepare the 
thermal state and probe to measure the residual rotation. 
Displace and measure (DM [a]): prepare the thermal state, 
displace by a and probe. Thanks to atom loss at each ther- 
malization, the atom number is varied by repeating DM sev- 
eral times. Measure number of atoms Na'- by pumping the 
atoms into F=l, mp=l and probing we measure the number 
of atoms in the trap. To correct for drifts, a sequence without 
displacement (DM[0]) is performed every 11 runs. We perform 
the sequence varying the displacement to acquire a dataset of 
quantum-noise-limited measurements of Pi NG '(Sy° ut ' ) ) for dif- 
ferent a. 



atoms, Sy 011 ^ is detected with a shot noise limited (SNL) 
balanced polarimeter in the ±45° basis. N-^ is measured 
with a beam-splitter and reference detector before the 
atoms. The probing-plus-detection system is shot-noise- 
limited above 3 x 10 5 photons/pulse. Previous work with 
this system has demonstrated QND measurement of the 
collective spin F z with an uncertainty of ~ 500 spins 
2EI2I. 



We generate Gaussian and non-Gaussian distributions 
with the following strategy: we prepare a "thermal state" 
(TS), an equal mixture of the F = \,mp = —1,0,1 
ground states, by repeated unpolarizcd optical pumping 
between the F = 1 and F = 2 hyperfine levels, finish- 
ing in F = 1 ,26|. By the central limit theorem, the 
TS of 10 6 atoms is nearly Gaussian with (F z ) = and 
var(F z ) = a 1 = 2^,4/3. By optical pumping with pulses 
of circularly-polarized light we displace this to (F z ) = 
a, with negligible change in var(i^ z ) j2?|, to produce 
P a (F z ) = {aV^)- 1 exp[-(i^ - a) 2 /(2a 2 )}. By displac- 
ing different TS alternately to a+ and a_, we produce 
an equal statistical mixture of the two displaced states, 
PL NG) {F Z ) = [P a+ (F Z ) + P a _(F z )}/2. With properly- 

chosen a± , Pq NG '' (F z ) closely approximates marginal 
distributions of mixtures of n = 0,1 Fock states and 
m . = N, N — 1 symmetric Dicke states. The experimental 
sequence is shown in Fig. [2] 

Detection, Analysis and Results - For each preparation, 
100 measurements of F z are made, with readings (i.e., 
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FIG. 3. (color online) Measured and predicted fc 4 with resid- 
uals for non- Gaussian distributions of different a. Readout 
noise is varied by the choice of Nr. Data is normalized to Nr 
and OA- Top: Points show normalized — k± calculated from 
N = 100 preparations of the ensemble with different a (hor- 
izontal axis), and Nr (colors). Black line indicates expected 
— «4, red line (top) shows — m + ^/var^) calculated from the 
distribution parameters for the largest readout noise. Bottom: 
normalized residuals ( — fc4 + k±) / y/ vax{k/±) . The normaliza- 
tion is done with the expected var(&4) for each Nr. Measured 
ki agrees well with theory, in particular, measurement noise 
increases the observed variance, but not the expectation. 

estimated F z values by numerical integration of the mea- 
sured signal) m, = 2S^ ut '^ /N'f 1 . Because the mea- 
surement is non-destructive and shot noise limited, we 
can combine Nr readings in a higher-sensitivity meta- 
pulse with reading M = X) TO i 1261 • This has the dis- 
tribution P a± (M) = cxp{-(M-a ± ) 2 /(2<7 2 M )]/(<7 M V2^) 
where the variance a 2 M = a\ N' A 2 N R -\-a 2 R includes atomic 
noise a 2 A N' A 2 and readout noise, a 2 R = Nr/Nl with 
N' A = N A /N A IAX . The variance a\ is determined from 
the scaling of var(M) with N a and Nr, as in [2(|. The 
readout noise can be varied over two orders of magnitude 
by appropriate choice of Nr. For one probe pulse and 
the maximum number of atoms we have ctr/ct\ = 84.7. 

To produce a non-Gaussian distribution, we compose 
metapulses from Nr samples drawn from displaced ther- 
mal state (DM[a + ] or DM[a_]) preparations with equal 
probability, giving distribution Pi NG) (M) = [P a+ {M) + 
P a _(M)]/2. With a M = (a+ - a_)/2, the distribu- 
tion has K2n+1 — 0, K2 = a 2 j + <j\j, k 4 — —2a\j, 
kq = 16a 6 M , Kg — ~212o? M . Our ability to measure the 
non-Gaussianity is determined by (/c 4 ) = k 4 and from Eq 
© 

var(fc 4 ) = 136Na s M /N {1) - U4N 2 a 4 M (a 2 M + a 2 M ) 2 /N {2) 
+24N 2 (N + l)(a 2 I + a 2 I ) 4 /N {3) . (5) 

As shown in Fig. [3l the experimentally obtained values 
agree well with theory, and confirm the independence 
from measurement noise. 



FIG. 4. (color online) Signal-to-noise in estimation of k 4 ver- 
sus readout noise for different a' = «m /{Nr a a). Points show 
measurement results, lines show theory, (details in the text) 

The "signal-to-noise ratio" for k 4 , S = K 4 2 /var(fc 4 ), is 
computed using Eq. ([5]). ka = —2a% r , and experimental 
otMi Nr, or, is shown as curves in Fig 4. We can confirm 
this <S* experimentally by computing Sjy = (fei) 2 /var(fc 4 ) 
using fc 4 values derived from several realizations of the 
experiment, each sampling P^ G N times. In the limit of 
many realizations Sn — > S. We employ a bootstrapping 
technique: From 100 samples of P^ G \m) for given pa- 
rameters aM, Nr and Na, we derive thirty-three TV = 20 
realizations by random sampling without replacement, 
and compute (k^) and var(fc 4 ) on the realizations. As 
shown in Fig. 01 good agreement with theory is observed. 
Optimum estimation of non- Gaussian distributions - Fi- 
nally, we note that in scenarios where measurements are 
expensive relative to state preparation (as might be the 
case for QND measurements of optical fields or for test- 
ing the successful storage of a single photon in a quantum 
memory), optimal use of measurement resources (e.g. 
measurement time) avoids both too few preparations and 
too few probings. 

We consider a scenario of practical interest for quan- 
tum networking: a heralded single-photon state is pro- 
duced and stored in an atomic ensemble quantum mem- 
ory. Assuming the ensemble is initially polarized in 
the X direction, the storage process maps the quadra- 
ture components X, P onto the corresponding atomic 
spin operators Xa,Pa oc F z ,—F y , respectively QND 
measurements of F z are used to estimate Xa, and 
thus the non-Gaussianity of the stored single photon. 
Due to imperfect storage, this will have the distribu- 
tion of a mixture of n = and n = 1 Fock states: 
p = (1 - p)|0)(0| + p|l)(l|. For a quadrature X, 
we have the following probability distribution P p (X) = 
cx.p[-x 2 / (2a 2 )] [px 2 j a 2 + 1 — p) j (\/27rao), where er is 
the width of the n = state. 

Taking in account the readout noise a\, the cumulants 
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FIG. 5. (color online) Signal-to noise-ratio S versus Nr for a 
fixed probe number NmNr = 1 x I0 5 for the probability dis- 
tribution associated with Fock state mixture described in the 
text with a normalized n=0 width <7o = 1. Red curve (top): 
p=l. Green curve (bottom): p= 0.5 with SNL measurement: 
(JR = ^20 /Nr. 



are K odd = 0, n 2 = (2p + l)a% + a%, k 4 = -I2p 2 af 1 , 
k 6 = 240p 3 crf f , k$ = -10080p 4 cr^, where the readout 
noise a R is included as above. Here K4 is directly related 
to the classicality of the state, since p > 0.5 implies a 
negative Wigncr distribution [19[. 

For a fixed total number of measurement resources 
NmNr, an optimal distribution of resources per measure- 
ment Nr exists as shown in Fig. [5j With increasing Nr, 
the signal-to-noise first increases due to the improvement 
of the measurement precision. Then, once the increased 
measurement precision no longer gives extra information 
about k&, the precision decreases due to reduced statis- 
tics because of the limited total number of probes. For 
a large total number of measurements, we can derive a 
simplified expression of this optimum. We derive asymp- 
totic expressions of S: Sl (Sh) for or <C <Jq (<jr 3> 00). 
The optimal Nr is found by solving Sl = Sr giving 
vr 8 ~ (?q( 1 + 8 P ~ 12 P 2 + 48 P 3 - 24p 4 ). For this optimal 
<jr, the measurement noise is in the same order of mag- 
nitude as the characteristic width of the non-Gaussian 
distribution. 

Conclusion - The cumulant-based methods described 
here should be very attractive for experiments with non- 
Gaussian states of material systems such as atomic en- 
sembles and nano-resonators, for which the state prepa- 
ration time is intrinsically longer, and for which measure- 
ment noise is a greater challenge than in optical systems. 
Cumulant-based estimation is simultaneously efficient, 
requiring few preparations and measurements, accommo- 
dates measurement noise in a natural way, and facilitates 
statistically-meaningful tests, e.g., of non-classicality. 
Experimental tests with a cold atomic ensemble demon- 
strate the method in a system highly suitable for quan- 
tum networking, while the theory applies equally to other 
quantum systems of current interest. 
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